Spin-state Crossover Model for the Magnetism of Iron Pnictides 



<N 

O 

<N 



& 
o 



> 



00 

o 

(N 



X: 



Jiri Chaloupka^'^ and Giniyat Khaliullin^ 

^ Max Planck Institute for Solid State Research, Heisenbergstrasse 1, D-70569 Stuttgart, Germany 
^Central European Institute of Technology, Masaryk University, Kotldfskd 2, 61137 Brno, Czech Republic 

(Dated: August 7, 2012) 

We introduce a minimal model describing magnetic behavior of Fe-based superconductors. The 
key ingredient of the model is a dynamical mixing of quasi- degenerate spin states of Fe^^ ion by 
intersite electron hoppings, resulting in an effective spin 5'eff in the ground state. The moments 
SeE tend to form singlet pairs, and may condense into a spin nematic phase due to the emergent 
biquadratic exchange couplings. We show that while the spin length 5'eff is robust against the 
variations of physical parameters, its long-range ordered part may take any value, resolving the 
puzzle of large but fluctuating Fe- moments observed. Underlying singlet correlations explain also 
the unusual temperature dependence of the paramagnetic spin susceptibility. 

PACS numbers: 75.10.Jm, 74.70.Xa, 71.27.+a 



Since the discovery of superconductivity (SC) in 
LaFeAsOi-icFa^ [1], a large number of Fe-based SC's have 
been found and studied in a great detail [2]. Evidence is 
mounting that quantum magnetism is an essential part 
of the physics of Fe-based SC's; in this regard, they are 
similar to heavy fermion and cuprate SC's. However, the 
origin of magnetic moments and the mechanisms that 
suppress their long-range order (LRO) in favor of SC are 
apparently different from Kondo or Mott physics that 
operate in rare-earth and cuprate compounds. 

The magnetic behavior of Fe-based SC's is unusual. 
The ordered moments range from 0.1 — 0.4 /iB, as in 
spin-density wave (SDW) metals like Cr, to 1 — 2 /iB 
typical for Mott insulators, causing debates whether the 
spin-Heisenberg or fermionic-SDW pictures [lol - 

Il4j are more adequate. The external/chemical pressure 
strongly affects the ordered moment values; however, 
irrespective to the strength or very presence of LRO, 
the Fe-ions possess universally the fluctuating moments 
^ 1 — 2 /iB jlll, [13] 5 even in apparently "nonmagnetic" 
LiFeAs and FeSe. In fact, ab-initio calculations suggested 
early on that the Fe- moments, "formed independently on 
fermiology" [l7^ and "present all the time" [3], are in- 
strumental to reproduce the measured bond-lengths and 
phonon spectra [1, [l7l - [l9| . Neutron scattering experi- 
ments [20] observe intense high-energy spin-waves that 
are almost independent of doping/ temperature, consis- 
tent with the picture of local moments induced by Hund's 
coupling [21] and coexisting with metallic bands. 

While the formation of the local moments in orbitally 
degenerate system is natural, it is surprising that these 
moments (residing on a simple square lattice) may or- 
der or remain disordered, depending on pressure, isova- 
lent substitutions, etc. Moreover, the Fe-pnictides are 
semimetals where the electron-hole pairs tend to con- 
dense into SDW state, further supporting magnetic order 
of the underlying moments. A fragile nature of the mag- 
netic order in Fe-pnictides thus implies the presence of 
a strong quantum disorder of local moments, not cap- 



tured in ab-initio calculations that invariably lead to 
large LRO-moments over an entire phase diagram. The 
ideas of domain wall motion [18] and local spin fluctu- 
ations [21] were proposed as a source of spin disorder, 
but no clear and tractable model of quantum magnetism 
in Fe-based SC's has emerged to date. Here we propose 
such a model. 

Since Fe-pnictides are distinct among the other (Mn, 
Co, Ni-based) pnictide families, their unique physics 
should be rooted in the specific features of the Fe-ion it- 
self. In fact, Fe^+ is famous for its spin-crossover [22*]: it 
may adopt either of 5=0, 1, 2 states depending on orbital 
splitting, covalency, and Hund's coupling. As the ionic 
radius of Fe is sensitive to its spin, Fe-X bond length 
{X denotes a ligand) is crucial and pressure reduces the 
spin value. In oxides, S=2 is typical and 6'=0, 1 occur 
at high pressures only (e.g., in the Earth's lower man- 
tle [23]). In compounds with more covalent Fe-X bonds 
(e.g., X=S, As, Se, . . . ), 5'=0 state is more common while 
6'=1, 2 levels are higher in energy. Here it comes the ba- 
sic idea of this Letter: when the covalency and Hund's 
coupling effects compete, the many-body ground state 
(GS) is a coherent superposition of different spin-states 
intermixed by electron hoppings, resulting in an aver- 
age effective spin 6'eff whose length depends on pressure, 
doping, etc. We design and solve a model exploring this 
dynamical spin-crossover idea, and find that: {i) inter- 
actions between S'eff contain large biquadratic exchange 
and favor singlet pairs, explaining unusual increase of 
the magnetic susceptibility with temperature 0, {ii) 
spin-nematic correlations emerge competing with mag- 
netic LRO, {Hi) the ordered moments m vary widely, 
< m < ^eff, but magnon spectra are universal and 
scale with Sq^ as in the experiment [lO, HH . We predict 
new collective (spin-length fluctuation) modes accessible 
by resonant x-ray scattering. 

The Fe-ions in pnictides have a formal valence state 
Fe^+((i^). Among its possible spin states [see Fig. [TJa)], 
6'=0 must have the lowest energy; otherwise, the ordered 
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moment would be too large and robust. The 15=0,1 
states, "zoomed-in" further in Fig. [Hb), are most im- 
portant since they can overlap in the many-body GS 
by an exchange of just two electrons between two ions, 
see Fig. [U^c). The corresponding a^— process converts 
Fe(5'=0)-Fe(5'=0) pair into Fe(6'=l)-Fe(5'=l) singlet 
pair and vice versa; this requires the interorbital hopping 
which is perfectly allowed for ~ 109° Fe-As-Fe bonding. 
Basically, hzis a. part of usual exchange process when local 
Hilbert space includes different spin states 5'=0,1; hence 
n ^ J . Coupling J between S=l triplets is contributed 
also by their indirect interaction via the electron-hole 
Stoner continuum. This contribution depends on the 
Fermi-surface topology and, as the band structure calcu- 
lations show [26] , reduces upon doping since the electron- 
hole balance of a parent semimetal is no longer perfect. 

The Hamiltonian describing the above physics com- 
prises three terms: on-site energy Et of S=l triplet T 
relative to 5'=0 singlet 5, and the bond interactions J: 



n- 



(ij) 



(1) 

The operator Djj creates a singlet pair of spinfull T- 
particles on bond (ij). For a general spin S of T-particles, 

D^J = Y.M{-^)^^^T,^+MTj^-M with M = -5, . . . , +5 

denoting the N = 2S ^ 1 projections; physically, N = Z. 
(The normalization factor l/\fN is left out for conve- 
nience). The constraint Usi + uti = 1 is implied. 

The above model rests on three basic features of 
Fe-pnictides/chalcogenides: {%) spin-state flexibility of 
Fe^+ that can be tuned by pressure increasing Et^ {H) 
edge-sharing FeX4 tetrahedral structure allowing "spin- 
mixing" /^-term, and (iii) semimetallic nature which 
makes J values to decrease upon doping [26|, [13 • 

Figure [TJd-f) demonstrates the behavior of the N = 3 
model (spin-1 T-particles) on a single bond. The GS 
wavefunction |V^gs) = cosq^IV^a) -\- sin a \iIjb) is a su- 
perposition of two singlet states \iIja) = ^jsj |vac) and 
\iPb) = :^^i2l^^^) iiiixed by the A^-term. The mixing 
angle is given by tan 2a = >/3/^/(£^t — J) and the GS en- 
ergy Egs = Et-J- ^{Et - J)2 + 3/^2. At /i: = 0, there 
is a sudden jump [Fig. Hl^e)] from tit = 0, 5* = state to 
TiT = 1, = 1 once the exchange energy compensates the 
cost of having two T-particles. The dynamical mixing of 
spin states due to n-ieiia converts this transition into 
a spin-crossover: the effective spin length 6'eff = ut in 
the GS increases gradually. Fig. [TJf) shows that /^-term 
strongly stabilizes the singlet pair of T-particles; we will 
see that this translates into a large biquadratic coupling 
(X [Si • which is essential in Fe-pnictides (26l. [joj. 

Turning to the model ([T]) on a square lattice, we notice 
first that for A/" ^00 and large the GS is dominated by 
tightly bound singlet dimers derived from the single-bond 
solution. The resonance of dimers on square plaquettes 
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FIG. 1. (a) Schematic view of low [S = 0), intermediate 
(5' = 1), and high (5' = 2) spin states of Fe^+(3(i^). (b) 5' = 
and 5* = 1 states differ in two electrons (out of six) occupying 
either the same or two different t2g orbit als. The S = 1 
state has a larger ionic radius, (c) The At-exchange process 
generating a singlet pair of 5* = 1 triplets T of two Fe^^ ions, 
both originally in the 5* = state (denoted by s). (d) The GS 
wavefunction of a Fe^^-Fe^^ pair is a coherent superposition 
of two possible total-singlet states, optimizing energy gain 
of the At- pro cesses, (e) Effective spin (average occupation of 
S = 1 state per Fe-ion), depending on the ratio of the coupling 
J between S — 1 states and their energy Et- (f) Energy levels 
labeled by the total spin value of the Fe^^-Fe^^ pair. Only 
singlet pairs are affected by With increasing At, the S = 1 
states are gradually mixed into the GS. 



then supports a columnar state [30] breaking lattice sym- 
metry without magnetic LRO. In the opposite limit of 
A" = 1, the model shows a condensation of hardcore T- 
bosons as n increases (and can be investigated using spin- 
wave approach [31]). We found that the A' = 3 model 
relevant here is also unstable (at sufficiently large J) 
towards a condensation of T-particles with S = 1. This 
condensate hosts interesting correlations not present in a 
conventional Heisenberg model. We discuss them based 
on the following variational wavefunction which describes 
Gutzwiller-projected condensate of spin-1 T-bosons: 

I*) = n [V^P 4 + (E dc^Tc^)^] I vac) , (2) 
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FIG. 2. (color online), (a) Condensate density p (= 5'eff) 
obtained from Eq. (|2]) as a function of angles i^, ip which 
parametrize the model ([T]) via Et = cosi^, ki = sin-z^cos 
and Ji — sini^ sin We set k^Ikx — J2/J1 — 0.7. (b) The 
ordered spin moment value m. (c) T-occupation per site tit 
obtained by an exact diagonalization, to be compared with p 
of panel (a). A bipartite 12-site cluster defined by the vec- 
tors (2,2) and (—4,2) in the square lattice was used, (d) 
Phase diagram and the ordered moment m as a function of Et 
and relative J-strength for fixed = 100 meV, Hi2 = O.Ta^i, 
= 140 meV, = 0.7j}°\ (e,f) Effective spin-length 
p — Sef£ and ordered moment m at the (e) Et = 100 meV 
and (f) J / J^^^ — 0.75 lines through the phase diagram in (d). 



where p G [0, 1] is the condensate density to be un- 
derstood as the effective spin length ^eff. The com- 



plex unit vectors di = Ui 



1) deter- 



mine the spin structure of the condensate in terms of 
the coherent states of spin-1 (sH, [33| corresponding to 
the operators Tq, {a = x^y^z): = (T+i — T_i)/>/2z, 
Ty = (T+i + T_i)/V^, = iTo. This is advantageous 
due to the symmetric expressions Dij = Ti^Tja and 
= —iea/s-fTpT^. The GS phase diagram obtained by 
minimizing and cross-checked by an exact diag- 

onalization on a small cluster is presented in Fig. [2j We 
have included nearest-neighbor (NN) and next-NN inter- 
actions and fixed their ratio at J2/J1 = 1^2/1^1 = 0.7, 
reflecting large next-NN overlap via As ions. Like in 
Ji — J2 model, this ratio decides between (tt, tt) and (tt, 0) 
order. Fig. [2ja,b) contains, apart from a disordered (un- 



condensed) phase {p = 0) at small hc^ J, three distinct 
phases depending on k/Et and J / Et values: {%) Ferro- 
quadrupolar (FQ) phase with Ui = u oi unit length and 
zero Vi. This phase is characterized by the quadrupo- 
lar order parameter {S'^S^ — ^S'^5a^) = p {^^ap — UaUs) 
with u playing the role of the quadrupolar director [33| , 
but it has zero magnetization. This state, often referred 



to as spin-nematic^ appears in biquadratic-exchange [32- 
im and spin-1 optical lattice models 36l-[39|. {ii) Non- 



saturated antiferromagnetic (ns-AF) phase with stripy 
magnetic order, specified by Ui = (0, 0, u) and Vi = 
(0,v,0)e*^-^^ with Q = (7r,0). The LRO-moment is 
given by m = 2puv which can take values from to 
'S'eff = p, even on a classical level, (iii) Saturated an- 
tiferromagnet (AF) with the same Q vector, but now 
with u = V = I/V2 and the ordered spin moment 
m = 5'eff = 1. 

The part of the phase diagram relevant to pnictides is 
shown in Fig. [2fd). The decrease of J is associated with 
doping that changes the nesting conditions [26.] , while the 
increase of Et is related to external/chemical pressure. 
Fig. [2fe,f) shows that the LRO-moment m quickly van- 
ishes as J {Et) values decrease (increase); however, the 
spin-length ^eff = p remains almost constant (~ 1/2), 
corresponding to a fluctuating magnetic moment ^ 1 pb- 
This quantum state is driven by /^-process which gener- 
ates the spin-1 states in a form of singlet pairs. 

We consider now the excitation spectrum, focusing on 
a realistic case of large condensate density (p ^ 0.4). It is 
convenient to separate fast (density) and slow (spin) fluc- 
tuations. To this end we introduce pseudospin r = 1/2 
indicating the presence of a T-particle, and a normal- 
ized vector field d defining the spin-1 operator as S' = 
—i{d^ X d) [classical part of d enters Eq. ([2])]. The Hamil- 
tonian then reads as 

^ = iSx ^ (i - rf ) - ^ {T+T+ di ■ dj + h.c.) 

i {ij) 

- E i\ - - (^I X d^) • id] X d,) , (3) 

(ij) 

and is decoupled on a mean-field level. The condensate 
spin dynamics is given by 0(3)-symmetric Hamiltonian 

= - E {di-dj^h.c)-Y^ Jij (4 X di ) • {d] x dj ) (4) 

with the renormalized kij = i^ij{T^Tj') ~ i^iji^ ~ p)p 
and Jij ~ Jijp^' The excitations above the GS ([2|) 
are found by introducing a, 6, c bosons according to 
d = {dx, dy, dz) = {a^ub — ive^^'-^c^ —ive^^'^h + uc)^ 
and replacing the condensed one as c, ^/T^^riZ^^rih. 
The resulting quadratic part of the a, 6 Hamiltonian is 
solved by the Bogoliubov transformation. A similar ap- 
proach is used for the r-sector Hamiltonian describing 
the condensate density fluctuations 5p. 
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FIG. 3. (color online), (a) Dispersion of the condensate den- 
sity (^p, solid-black) and the ordered moment- length (^|m|, 
dotted-blue) fluctuations, and the magnon dispersion (solid- 
blue), at the point A in the phase diagram of Fig. [2jd). All 
3 modes are active with respect to resonant x-ray scatter- 
ing, and the latter 2 to neutron scattering, (b) Evolution of 
the magnetic excitations going from FQ to the ns-AF phase 
[C B ^ A in Fig. [2{d)]. Two-fold degenerate quadrupole- 
waves (C) split into the magnon (solid lines) and the ^|m| 
mode (dotted lines). The latter represents oscillations be- 
tween the nematic and magnetic orderings and is gapful. 



Shown in Fig. [3] is the dispersion of the excitations 
for several points of the phase diagram. The density 
(i.e., ^eff) fluctuations are high in energy. Fig. [3jb) 
focuses on the magnetic excitations. In the FQ phase, 
quadrupole/magnetic modes are degenerate and gapless 
at gr = 0, where they correspond to the Goldstone modes 
associated with a free director rotation. As the ns-AF 
phase is approached, the gap at Q decreases, and closes 
upon entering the magnetic phase. Importantly, the spin 
fluctuation spectra is determined by the effective spin 
5'eff = p and not by the ordered moment m value. Since 
p varies only slightly, the magnon energies should be com- 
mon to different materials, as in fact observed 20|, 25| . 

The magnetic modes in Fig. ^h) resemble excitations 
of bilinear-biquadratic spin model [33]. In fact, the dis- 
persion in FQ phase can be exactly reproduced [40] from 
an effective spin-1 model XI ^^j) Jij^i • — hij{Si • Sj)'^^ 

with J and k shown above. A large biquadratic coupling 
was indeed found to account for many observations in 
iron pnictides d, [26|, [Soj. We note however, that this 
model possesses FQ and AF phases only and misses the 
ns-AF phase, where the ordered moment is reduced al- 
ready at the classical level; also, it does not contain the 
key notion of the original spin-crossover model, i.e., for- 
mation of the effective spin S'eff and its fluctuations. 
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FIG. 4. (a) Density of states of the magnetic excitations cal- 
culated for the point C of Fig. [2jd). We included the damp- 
ing (e.g., due to coupling to the Stoner continuum) in a form 
r(cj) min(a;,r) with F ujq/2. The result with F = is 
shown for comparison, (b) Temperature dependence of the 
uniform susceptibility x- The components x\\ (x±) parallel 
(perpendicular) to the local director u are also shown. Low- 
energy cutoff of 1 meV was used. 



Singlet correlations inherent to the model may ex- 
plain also unusual increase of the paramagnetic suscep- 
tibility x(^) with temperature [24]. Considering non- 
magnetic FQ phase, we find that for the field parallel 
to the director n, x is temperature dependent, X|| = 
2^ / dcjA/'(cj) sinh"^ ^, where J\f{uj) = ^(^ - ^q) 
is the density of states (DOS) of magnetic excitations, 
while x± for the field perpendicular to u is constant and 
inversely proportional to the bandwidth of excitations. 
(The physical x contains additional factor of p^g'^p^NA 
with g = 2 and the Avogadro number Na)- Their average 
X = (X|| +2x±)/3 with respect to the local director orien- 
tation corresponds to the measured x{T)^ assuming slow 
rotations of the director. The DOS shown in Fig. Hl^a) 
is contributed mainly by the regions around (tt, 0) and 
(0,7r) where AF correlations do reside. The correspond- 
ing thermal excitations lead to the increase of x 
very high temperatures [see Fig. HJ^b)]. 

To conclude, we proposed the model describing quan- 
tum magnetism of Fe-pnictides. Their universal magnetic 
spectra, wide-range variations of the LRO-moments, 
emergent biquadratic-spin couplings are explained. The 
model stands also on its own: extending the Heisenberg 
models to the case of "mixed-spin" ions, it represents 
novel many-body problem explored here only in part and 
deserves future study. Of a particular interest is the effect 
of band fermions (only mentioned above as the origin of 
doping dependent J values and magnon damping) which 
should have a strong impact on low energy dynamics of 
the model, e.g., converting the gr = Goldstone modes 
into overdamped spin- nematic fluctuations. Understand- 
ing the effects of coupling between local moments and 
band fermions, including implications for SC, should be 
the next step towards a complete theory of Fe-pnictides. 
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